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ABSTRACT 

Radio observations of Fanaroff-Riley class II sources often show correlations between 
the synchrotron emission and the linear-polarimetric distributions. Magnetic position 
vectors seem to align with the projected emission of both the radio jets and the sources' 
edges. Using statistics we study such relation as well as its unknown time evolution 
via synthetic polarisation maps of model FR II sources formed in 3D-MIID numerical 
simulations of bipolar, hypersonic and weakly magnetised jets. The magnetic field is 
initially random with a Kolmogorov power spectrum, everywhere. We investigate the 
structure and evolution of magnetic fields in the sources as a function of the power of 
jets and the observational viewing angle. Our synthetic polarisation maps agree with 
observations, showing B-field vectors which are predominantly aligned with the jet 
axis, and show that magnetic fields inside sources are shaped by the jets' backfiow. 
Polarimetry is found to correlate with time, the viewing angle and the jet-to-ambient 
density contrast. The magnetic structure inside thin elongated sources is more uniform 
than inside more spherical ones. We see jets increase the magnetic energy in cocoons 
in proportion to the jet velocity and the cocoon width. Filaments in the synthetic 
emission maps suggest turbulence develops in evolved sources. 

Key words: galaxies: jets - galaxies: active - intergalactic medium - methods: 
numerical - MHD - turbulence 



1 INTRODUCTION 

Centimetric wavelength observations reveal synchrotron 
emission from extraRalactic F anaroff-Riley Class II radio 
sources (iFanaroff fc Rilevll 19741; FR lis hereafter) and radio- 
loud quasars (see e.g. iBridle fc Perle v 1984 and references 
therein). Linear polarisation fractions within ~ 10-50 % 
are commonly seen in these objects. Polarisation maps 
of these sources show patchy distributio ns which corre- 
late with the luminosity distribution (see ISaikia fc Salter! 
Il988l for a review). Projected magnetic field vectors 
are predominantly parallel to both the radio jets and 
the bound aries of radio lobes (Alexander, Brown & 
Scott 1984: iBridle fc PerlevI Il984l : Leahy, Pooley & Riley 
1986, 1997; 'Black et al."l992; Johnson, Leahy fc Garring- 
ton 1995; iHardcastle et al.. ■1997. . 1998; iLeahv et all 1 19971 ; 
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[Gilbert et al]|2004l : lMullin et al.ll2006l ). Strong emission gra- 
dients are often followed by the vectors perpendicularly, and 
when multiple hotspots are observed in one of the two radio 
lobes, the vectors seem to follow a line that would connect 
the hotspots. Linear polarisation fractions of radio jets tend 
to be higher at the edges than at inner regions. Linear po- 
larisation fractions of radio lobes tend to be higher at the 
edg es than at regions both inside and between the lobes 
(tSaikia fc Salteijll988l ). 

The direction of the magnetic field component that is 
in the plane of the sky is often inferred by computing the 
Stokes parameters on the observed signal. It is possible to do 
the calculations inversely in order to model the polarimetry 
distri bution that r esults from given magnetic field geome- 
tries (iLaind Il98ll ; iJoned Il988l ) . Such studies indicate that 
magnetic fields in FR lis seem to consist of a combination 
of ordered and disordered (anisotropic) fields along the jets 
and their vicinities, as wel l as a random component at the in- 
ner regions of radio lobes l|Laindll98ll : IBridle fc Perievlll984l ; 
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ISaikia fc Salted 1 19881 ). Circumferential magnetic structures 
are also frequently observed in the outer ede e of the sources 
ijBridle fc Perlevl(l984l : ISaikia fc Saite3ll988h . 

Based on the radio lu minosity distribution observed in 
several FR IIs ,'Recs' ('l97l|), Lo ngair, Ryle fc S cheuer (1973), 
iBlandford fc^Rees (, igT^ and [Scheueif (| 19741 ) proposed the 
following model for the plasma dynamics in the sources. 
Magnetised relativistic plasma jets are launched from a cen- 
tral engine located inside active galactic nuclei (AGN) which 
are typically seen at positions that match those of the ra- 
dio cores. A cavity (the cocoon, hereafter) is inflated with 
the jets' plasma and a strong bow shock is driven on the 
ICM. Jets collide with the ambient medium at their work- 
ing surfaces. Radio hotspots are seen at the leading edges of 
the lobes because the plasma pressure is the highest there. 
The plasma nearby is pushed towards the radio core and a 
backflow of magnetised plasma develops. Radio synchrotron 
lobes are thus formed, separated from the ICM by a contact 
discontinuity. 

At kiloparsec-scales, magnetic fields in FR lis are of- 
ten modeled in energy equipartition with the synchrotron 
emitting ultra-relativistic electrons. Magnetic flux freezing 
is expected to bond field lines and radio emitting electrons 
dynamically, hence the jets' backfiow should play an i mpor- 
tant r ole in shaping magnetic fie l ds in si de cocoons |Laing 
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[al 

1984 [Mnier 1985; Saik ia fc Saiteii r 

The expansion of radio sources must be considered 
to understand their inferred magnetic structure. Evolu- 
tionary models have provided analytical expressions for 
the global time dependence of the volume, the pressure 
and the energy inside cocoons (S cheuer 1974; Fallc 199ll: 
Begelman fc Cioffil IT989I : iKaiser fc Alexander! Il997l : iKrausd 
20031 : Heinz. Reynolds fc Begelman 1998). The large- 
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scale features of the complex non-linear dynamics in such 
plasma cavities have been captured by numerical simulations 
(for r eviews see | Norman 1993; Ferrari 1998; Pudritz ct ah, 
l2007h . Two-dimensional axisymmetrical simulations of mag- 
netised jets have confirmed the basic picture regarding jets 
(beams), lobes (cocoons) and bow shocks (e.g. Clarke, Nor- 
man fc Burns 1986; Kossl, Miiller fc Hill ebrandt 1990a, b ; 

1989: Frank et al. 1998; Stone fc Harded l2000l : 



Lind et al] 



Komissarov 19991 ) . These simulations also show that cocoons 



consist of a series of vortices. These structures arise in a com- 
plex feedback loop, where pressure modulations in the co- 
coons interact with the beams' shock pattern which, in turn, 
modifies the vortex shedding. The vo rtices decay in a tur- 
bulen t cascade (compare section 3 in iKrause. fc Alexander] 
l2007h . The latter process results in some degree of isotropi- 
sation of the field lines which should affect the alignment of 
magnetic fields and the fractional polarisation. 

Further, the expansion of cocoons involves magnetic 
field ampl ification. This happens via t wo field line stretching 
processes l|Matthews fc Scheueij[l990l ) : The poloidal stretch- 
ing mechanism, which arises because the fluid elements in 
the beam located close to the beam boundary take small 
turns, and thus end up towards the inner part of cocoons. 
In contrast, the fluid elements close to the jet axis make 
larger turns and end up near the outer cocoon boundary. Be- 
cause of the larger path length of these fluid elements, they 
lag behind. Hence shear amplifles the magnetic fleld in co- 
coons along the direction of the jet axis. On the other hand. 



the toroidal stretching process amplifles the toroidal compo- 
nent of magnetic flelds via cocoon expansion perpendicular 
to the jet axis. Unless the flow structure is axisymmetric, 
which is an unlikely configuration for real radio sources, the 
toroidal magnetic field may again be sheared, and thereby 
converted into poloidal field. To first order, the resulting 
magnetic field structure is determined by both, these com- 
peting processes and the initial condition. This picture has 
been refined by Gaibler, Krause fc Camenzind (2009) who 
initialised their simulation with a helical magnetic field con- 
fined to the beam. The poloidal component of these flelds 
returns to the source along, and close to, the beam, and 
therefore its strength drops steeply with distance from the 
jet axis, R. The radial cocoon expansion puts work into the 
toroidal fleld w hich, consequently, increases linearly with R, 
as predicted bv lMatthews fc Scheuerl l|l990l . toroidal stretch- 
ing). Hence, the magnetic energy in radio lobes could be 
largely due to dynamo action in the large scale jet flow, with 
little dependence on the conditions (set up) at the base of the 
beam. The literature on three-dimensional jet simulations, 
in contrast, has not paid much attention to these issues. 
In general, one flnds jet instabilities which are transparent 
to simulations with less degrees of freedom, e.g. je t fluting, 
deflec tion, disconnection and splash-back (see e.g. iNormar] 
Il993l ). Jet propagation has been studied also with relativistic 



MHD codes Ce.g. lLeismann et al.1120051 : iKeppens et~al]|2008l : 
iMignone et al.ll2010l . and references therein). These studies 
have not particularly focussed on polarisation properties. 
Relativistic jets have narrower cocoons for a given rest mass 
density ratio, and more stable beams. The motions in radio 
lobes are subrelativistic, and hence their physics should not 
be too much influenced by a relativistic nature of the jet. 

Synthetic observations are produced using data from 
numeri cal simulations in orde r to co mpare them with obser- 
vations. [MatthewT^fc^^Scheue^ l|l990l ) simulated the hydrody- 
namic advection and polarized synchrotron emission of ran- 
dom, pasive magnetic fields in AGN j ets, find i ng hig h linear 
polarisation fractions, of about 70 %. Iciarkel l| 19931 ) carried 
out 3D simulations of the interaction of a jet and a cloud 
with passive uniform magnetic fields. The synthetic syn- 
chrotron emission ma ps of Clarke sho wed filaments, formed 
by velocity shear. Hardee, C larke, fc Rosen (1997) carried 
out 3D-MHD simulations of supermagnetosonic magnetised 
perturbed equilibrium beams, where a section of an infinitely 
long beam is studied, and found synthetic intensity struc- 
tures similar to the ones observed in the jets of Cygnus A. 
More recently, Tregillis, Jones fc Ryu (2004b) investigated 
the fractional polarisation of synthetic synchrotron observa- 
tions of 3D-MHD AGN jet simulations. They found rather 
high fractional polarisations in regions where shock acceler- 
ation increases the emissivity, but much smaller fractional 
polarisation at regions where relativistic particles illuminate 
the volume more uniformly. In general, little attention has 
been given to the statistics of synthetic polarimetry and the 
way it relates with the properties of radio jets. 

In this paper we present 3D-MHD numerical simula- 
tions of hypersonic magnetised jets as well as synthetic 
synchrotron and polarisation observations. In contrast to 
Gaibler et al. (2009) and much other work, we do not start 
with a regular magnetic field component within the jet, but 
rely entirely on the field amplification due to the dynamics 
of cocoons (compare above) to create structure. Regarding 
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analysis and the ques t ions w e address, we follow essentially 
iMatthews fc Scheueij (|l990l ) with the important improve- 
ment that here we use a full three-dimensional magneto- 
hydrodynamic treatment for the jet simulation. 

This paper is organised as follows: in Section [2] we de- 
scribe the formalism of ideal MHD and the numerical meth- 
ods we use. Our implementation of the ICM, CMFs and 
AGN jets are also described there along with details of our 
calculations for the synthetic synchrotron emission and po- 
larimetry. In Section [3] we talk about the flow structure 
in our model sources and analyse it in terms of energet- 
ics. Synthetic maps are then presented and compared with 
FR II radio observations. The results are then interpreted 
and analysed statistically. Section |4] is dedicated to compare 
our models with previous numerical simulations. We then 
summarise and conclude our study in Section [5] which is 
followed by the bibliography. 



2 SIMULATIONS 

We describe the dynamics of plasma in the ICM and AGN 
radio jets using the system of nonlinear time-dependent hy- 
perbolic equations of ideal compressible MHD. In three di- 
mensions and non-dimensional conservative form, these are 
given by: 



dp 
dt 



+ V ■ (pV) 



+ V ■ (pVV -t- + 5^2 - BB) 



= ft (1) 
= pg + Pj(2) 



dE 



^+V - [(iJ + p, + BV2)V-B(V.B)] 



f -Vx(VxB) 



= 0, 



(3) 
(4) 



where p, Pg, V and B are the plasma density, thermal 
pressure, flow velocity and magnetic fields, respectively. In 
1(3)), S = pg/(7 - 1) -I- pV^/2 -I- _B^/2 and represents the to- 
tal energy density whereas 7 is the ratio of specific heats. 
In the right hand side of (O, Q and ([2} source terms are 
used to implement jets by injecting mass, pj, momentum, 
Pj , and kinetic energy, Ej (see Section 12. 2p , as well as a 
Newtonian gravitational acceleration, g, to keep the plasma 
in magneto- hydrostatic equilibrium (see Section [2. l.ip . 

We solve the above equati ons in three di mensions using 
the numerical code Flash 3.1 (|Frvxell et all [2000 ). Flash's 
new multidimensional unsplit constrained transport solver is 
employed to maintain th e divergence of magnetic fields down 
to <10~^^ (iLee fc Dean e 2008). A diffusive HLLC solver 
IIBatten et al.l 119971 ') prevents spurious low pressure and den- 
sity values from appearing in the grid. We use a Courant- 
Friedrichs-Lewy parameter of 0.25 and periodic boundary 
conditions in all the domain's faces. These boundary condi- 
tions prevent numerical noise from polluting the turbulent 
magnetic spectrum in the grid fSection l2.1.2l) . Our compu- 
tational domain is a cube with edges |x| 1/2, in compu- 
tational units, and has a uniform grid with 200'^ cells. This 
represents a volume of 200 kpc'^ meant to simulate the core 
of a cluster. 

We carried out five jet simulations (see Table 1) de- 
signed to experiment with the power of jets in terms of 



their velocities and densities. Computations ran for approxi- 
mately 12 hours on 64 processors at the CamGricj^ cluster of 
the University of Cambridge, and the production runs exe- 
cuted for about 4 hours (using 64 processors) at the Darwir0 
supercomputer of the University of Cambridge HPC facility. 



2.1 Initial conditions 

2.1.1 The ICM 

The cluster plasma is implemented using an equation of 
state of an ideal monoatomic gas, with a ratio of specific 
heats 7 = 5/3, a constant sound speed (c^ = "iVgl P = 1) 
thro ughout the domain and a density following a King pro- 
file (|King|ll972l ') 



PiCM(r) 



(l + (r/rc)2)3/J/2' 



(5) 



where the central density, pc, the central radius, r^, and /3 
take the values of 1, 0.8 and 2/3, respectively. 

To keep the magnetised gas in magneto-hydrostatic 
equilibrium we implement a radial acceleration source term 
g to equation and take the balance between this term 
and the total plasma pressure pg + /2. In the radial di- 
rection this term takes the form: 



9r 



2ci 



^rl (1 -I- (r/rc)2) 



1 + 



(6) 



where /3ni is the ratio of thermal pressure to magnetic pres- 
sure. 



2,. 1.2 Cluster magnetic field 

The magnetic field within the cluster is set up as an isotropic 
random field with a power law energy spectrum. Following 
iTribbld l|l99ll) and lMurgia et al.l (|2004h we generate a cubic 
grid in Fourier space, with 200'^ cells. For each of these, we 
define three components of a vector potential which takes 
the form A(k) = A(k)e**'''-', where k is the frequency vec- 
tor (fc^ — k'^ + ky + k^), i is the unitary complex number, 
while A and 9 are the vector amplitudes and phases, respec- 
tively. We draw 9(k) from a uniform random distribution 
within and 2tt, and A(k) is also randomly distributed but 
has a Rayleigh probability distribution 



P{A,e)dAde = 



A 



2TT\Ak\ 



-exp 



A' 



2\Ak\ 



dAde, 



where we choose the power law Ansatz 



(7) 



(8) 



for a given slope (. 

We transform to real space by taking the inverse fast 
Fourier Transform (Press ot al. 1992) of A(k). The resulting 
magnetic vector potential, A(x), is multiplied by the plasma 
density radial profile ([5]). This product implements mag- 
netic fiux freezing by generating fields, the strength of which 
follows the plasma density, and pressure, profile. The com- 
ponents of the vector potential are then read and mapped 



^ http:/ /www. cscience.cam.ac.uk/projects/camgrid/ 
2 http:/ /www. hpc.cam.ac.uk/darwin. html 
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into the staggered-grid cell interfaces of FlashS.l and the 
curl of this vector is then calculated to give the magnetic 
field. Finally, we normalise the resulting field so that the 
ICM's thermal pressure is approximately ten times larger 
than its magnetic pressure (/3m = Pg / {B'^ /2) ~ 10) every- 
where in the grid, whic h is a reasonable value in this context 
jCariUi fc Tavloij|2002l '). 

This procedure yields solenoidal magnetic fields tangled 
at scales of order our computational resolution and char- 
acterized by spatial variations following a magnetic power 
spectrum with a power law of the form 

(X = k-", (9) 

where we choose a Kolmogorov three- dimensional turbulent 
slope n = -11/3, based on the w ork of lVogt fc EnfilinI (|2003l . 
I2OO5I ) and lGuidetti" et alJ (l2008l ). We use the same realisation 
for all our runs. We note that the Fourier method implicitly 
imposes maximum and minimum scale on the field. 

We let this plasma relax for one crossing time before 
injecting the jets. 

2.2 Jets 

By implementing source terms to equations ([l|, ((2| and ((3]) 
we inject mass, a;— momentum and kinetic energy to the 
central grid cells that are within a control cylinder of ra- 
dius rj and height hj, resolved by 3 and 8 cells, respectively. 
Inside this "nozzle" we update the plasma density and x- 
velocity via constant source terms pj and Vj . Jets are contin- 
uously injected until they reach the computational bound- 
aries and then the simulations are stopped. Plasma pres- 
sure in the nozzle, pj, takes the constant value of the cen- 
tral ambient pressure (i.e. Pc/7). The jet density is com- 
puted using Pj = rjpc where the parameter rj takes the 
(low) values given in Table 1. We assume an ideal gas 
equation of state with 7 =5/3 for the jet material. The 
light densities of ou r jets are motivate d by t he work of 
lAlexander fc PoolevI (|l996l ') and iKrausd (|2003l ). and their 
high Mach numbers are based on the observed jets sided- 
ness associated with Doppler beaming, suggesting that FR II 
sources are at least close to relativistic up to scales of or- 
der 100 kpc ijMuUin fc Hardcas tle 2009 ). The Mach numbers 
of our jets with respect to the sound speed in the ambi- 
ent gas are 40, 80 and 130 which correspond to velocities 
close to 66x10^, 133x10^ and 216x10^ kms"\ respectively. 
The Mach numbers of the jets with respect to the sound 
speed in the beam material are 2.5, 5.7, 5, 11.3 and 8.2, 
as they appear i n Tab le 1. We extend the implementation 
of[0 mma et "aP (|2004h to simulate bipolar magnetised jets. 
The launch and collimation of the jets are assumed to oc- 
cur in the AGN "central engine" located at sub-resolution 
scales. 

The initial jet magnetic fields are kept from the ini- 
tialisation of the ambient medium, and no magnetic source 
term is applied. It therefore has a random topology, an aver- 
age /3m ~ 10 and, given the assumed power spectrum (Sec- 
tion [27lT2]), it is fairly uniform at scales ~ rj. We note, how- 
ever, there is no reason to believe the magnetic fields in FR II 
radio jets are related to the CMFs near the AGN; jet fields 
are expected to be advected up the beam from the central 
engine. Our choice of initial jet magnetic fields is based on 
the fact they seem to be weak at kiloparsec-scales and to 



Table 1. Simulations and parameters. 



Simulation 










name 


[Macii] 




[xlO^* W] 


[Myr] 


lighter-slow 


40 


0.004 


4.6 


14.1 


light-slow 


40 


0.020 


17.2 


8.3 


lighter-fast 


80 


0.004 


28.1 


7.1 


light-fast 


80 


0.020 


128.8 


4.4 


lighter-faster 


130 


0.004 


112.8 


4.7 



Time average average jet velocity in the nozzle. It is equal 
to the external Mach number. 

^ Time average average jet to ambient density contrast in the 
nozzle. 

Jet power from equation l lllll . 
Simulation end times. 

have a random component (Section 1). This is the case of 
the central fields in our model. 

As our jets propagate, their magnetic fields are de- 
formed by shear. The time averaged average beam /3ni is 
of about 50. The power of jets is the sum of thermal and 
kinetic power terms: 




which takes the following form at the grid 

= ^ivPc){'^r])vf + pj(TTrj)vj. (11) 



2.2.1 Cocoon contact surface 

We use a passive incompressible tracer, r(x, t), which is in- 
jected with the jet plasma to distinguish it from that of the 
ambient medium. When jet injection starts, r(x, tjet =0) 
takes the values of 0.99 and 1x10"^" at the nozzle and at 
the ICM, respectively. The tracer is then advected with the 
jet gas and takes values within 1x10"^° to 0.99. A com- 
parison of the distributions of r and p allows us to identify 
the contact surface of the cocoon with an accuracy up to 4 
computational cells. 

2.3 Synthetic radio maps 

Our simulations produce three-dimensional data cubes with 
information about the distribution of the magnetised gas in 
our model sources at different times during their expansion. 
Synthetic synchrotron emission and polarimetry are com- 
puted under the assumption that the radiation is linearly 
polarized. Beaming and light-travel effects are assumed to 
be negligible. Synthetic observations are produced at view- 
ing angles, 9^, of 30, 60 and 90 degrees measured from the 
jet axis to the line of sight (thus jets are in the plane of 
the sky when 9^ =90°). Given a viewing angle and a sim- 
ulation timestep, t, Stokes parameters are (i) calculated for 
every computational cell inside the source, using the mag- 
netic field components in the plane of the sky, B'(x,f), (ii) 
integrated through the source, along the line of sight, (t). 
Mathematically, 
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Mach=80, r] =0.004 & tj^t =7.1Myr Mach=80, r] =0.02 & tj^t =4.4Myr 



tutu 
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aim 4daa tri« 



atm tttu tm iM* 
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-4,1 $ft 9.1 #.i 



Figure 1. Two dimensional cuts (z =0) of the density (top row), the pressure (middle row) and the magnetic field strength (bottom 
row) distributions. The Mach number and density of the runs are given on the top of each column. Colour scales are logarithmic and 
show variables in corresponding computational units. We see a clear relation between the intrinsic structure of sources and the resultant 
field structure in the synthetic polarisation maps (Figures [4}j8j . 
© 2009 RAS, MNRAS OOO.flti^TI 
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I fJA^ 

fh) 'thKa6tii\ aver Louil n^u^cLic ^id (teft} mid uxiuJ [j^tt toroiJrkl imigDL'^.ijc fxLd (Ti^hi). 17 =£i.0OI (tap) 
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Figure 2. Two dimensional cuts (2 =0) of the magnetic fields for the Mach 80 jets (see Table 1). Colour scales are linear and show 
variables in corresponding computational units. 



/(xx,t) = 5(r)r(x,f)p,(x,f)[Bi(x,t)2+ where 

Q(xx,t) = ^i^S/g S{T)r{^,t)p4x,t)[BU^,tf- ,^2) 

Ui^^,t) = 5(r)r(x,t)p,(x,f)2B;(x,t) ^ir) ^ 

B'y{-K,t)dZ\\{t), 



1 for r(x,t) G [0.5, .99] (cocoon); 



(13) 



for r(x,t) e [1 X 10"", 0.5) (ambient), 
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: 80,11=0.020 
: 80,Ti=0.004 
■■ 40,ri=0.020 
40, Ti=0.004 
:130, Ti=0.004 




0.8 1.0 1.2 1.4 
Source size [100 kpc] 

Figure 3. Time evolution of the magnetic energy in the cocoon. 
Only the 7] =0.004 sources show flat gradients, typical of MHD 
turbulent dynamos. 



and x^, pc(x, t) and / represent the coordinates in the plane 
of the sky, the distribution of the cocoon pressure and the 
total intensity of the radiation, respectively. We note (|12|) 
are valid for a synchrotron emission spectral index a = 1, 
yet the degre e of polarisa tion we predict does not vary too 
much with a fLaing" 1980). The factor of 0.75 in the expres- 
sion of Q and U in p2p accounts for the maximum degree of 
linear polarisation for a uniform magnetic field and a power- 
law electron energy distribution. We model the density dis- 
tribution of synchrotron emitting electrons via the factor 
t(x, t) pc(x, in (|12[) . We do not follow any explicit energy 
gain or loss processes such as synchrotron cooling or shock 
acceleration (i.e. the background plasma pressure is propor- 
tional to the constant factor in the energy distribution of 
relativistic electrons). A detailed treatment of the electron 
distribution is beyond the scope of this paper. We note that 
we have also tried a constant density of radiating electrons, 
which did not significantly change the results. The polari- 
sation angle of the magnetic vectors, xb, and the degree of 
linear polarisation {— fractional polarisation), p, are given 

by 



1 TV 

Xb = - arctan(J7/Q) + -, 



P = 



(14) 



3 RESULTS AND ANALYSIS 

Synthetic polarisation and emission maps are presented in 
pairs characterized by the jet velocity (same as the external 
Mach number), the density contrast rj, the viewing angle 9^ 
and the time that jets have been active, tjct. Polarisation 
maps have a constant vector density of 0.5cells~'^. 



3.1 Flow structure 

The hydrodynamic flow structure of our simulations is very 
similar to what is generally found in the literature (compare 
section[TJ above, Figure[T}. The hypersonic jets fiow straight 
for a certain distance (2D-phase). Then three-dimensional 
instabilities develop, more clearly in the lighter jet runs, and 
the jet head oscillates around the jet axis, consistently with 



Scheuer's dentist drill (3D-phase). Cocoons are wider for 
lower jet density and faster jets, as expected. The relatively 
heavier jets {ri — 2x 10~^) propagate faster in the axial direc- 
tion, and backflows in their cocoons are much less turbulent 
than in their relatively lighter (77 =4xl0~^) counterparts. 

The evolution of cocoon magnetic fields is driven by 
the following dynamics. The field is initially random inside 
the jet injection volume. The injected momentum stretches 
field lines along the jet direction. This puts energy into 
the axial field, which is therefore amplified. The other field 
components are simply advected out of the injection vol- 
ume, and their field strength drops with time, within the 
injection volume and the beam. This process results in a 
poloidally dominated magnetic field, s imilar to the setup in 
iGaibler. Krause. fc Camenzindl lj2009l ). yet with some im- 
portant difi^erences: First, the field in the axial direction is 
patchy (compare Figure [2] for this and other details of the 
magnetic field), i.e. adjacent parts of the beam have the 
field parallel and anti-parallel to the fiow direction. Second, 
for a given plane perpendicular to the fiow vector, the field 
may in principle also be patchy, i.e. there is not necessarily 
one dominant toroidal field loop, but possibly two or more 
field loops across a section of the jet. However, the fact that 
the power spectrum used for the initial field setup strongly 
favours larger scales, still produces a predominantly toroidal 
configuration for the magnetic field perpendicular to the jet 
axis. 

As Gaibler et al. (2009) do, we find that the axial field 
lines return to the injection region very close to the jet. In 
our case this may even happen inside the beam, since any 
beam cross section may in general contain axial field patches 
of opposing directions. In the presence of a backflow, this re- 
quires field line reconnection, which should be easily possible 
in the jet head on numerical grounds due the complex flow 
pattern in this region. This seems to suggest that the gen- 
eral structure of the poloidal magnetic fleld does at least 
not very much depend on the initial conditions. The reason 
is that the elongation of the beam stretches the axial field 
lines and therefore amplifies the axial field. To have the field 
lines going out in the beam and returning close to it or even 
within is the configuration which requires the least amount 
of energy, and is therefore chosen by the system. 

Gaibler et al. (2009) find the toroidal part of the field, 
which cannot be lost to other field components via turbu- 
lence due to the axisymmetry condition in their study, in- 
creases linearly with distance from the jet axis. This may be 
easily understood from the induction equation. The phys- 
ical reason is the work done by the expanding cocoon on 
the toroidal field component is stored in that part of the 
magnetic field. We do not observe such a linear increase in 
the toroidal field in our 3D simulations directly, but we ex- 
pect this process also to be at work. Since it is related to 
the cocoon expansion, we expect more magnetic energy to 
be created by fatter cocoons, i.e. for lower jet density (most 
important for the cocoon width), and higher jet velocity, 
which is also found by Gaibler et al. (2009). They also show 
that this process is able to enhance the total magnetic en- 
ergy in the jet by a factor of a few (see their Figure 20). 
For our simulations, therefore, we expect a noticeable in- 
crease in the magnetic energy during the simulation time; 
the fatter the cocoons the higher the energy rise. Figure [3] 
shows this expectation is exactly what we find: Here we plot 
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Figure 4. Synthetic observations of the source with rj =0.004 and Mach = 40. Left: polarisation maps. Vectors follow XB and their 
length is given by p. Vectors are superimposed on linear contours of // (/}. Right: logarithmic grayscale maps of // (/). 
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Figure 5. Synthetic observations of the source with ri =0.004 and Mach = 80. Left: polarisation maps. Vectors follow XB and their 
length is given by p. Vectors are superimposed on linear contours of // (/). Right: logarithmic grayscale maps of // (/). 
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Figure 6. Synthetic observations of the source with -q =0.004 and Mach = 130. Left: polarisation maps. Vectors follow XB and their 
length is given by p. Vectors are superimposed on linear contours of // (/}. Right: logarithmic grayscale maps of // (/). 



© 2009 RAS, MNRAS OOO.mfSTI 



3D-MHD simulations of evolving magnetic fields in FR II radio sources 11 




fn) Time cstiliitiun m tii,- = 90" . F^nm lup to baK tdm t.„r = {'3 .7, -^^.'i. H. [ -t 




I 75 ^, veclprg 



0.50, i.CK). 7.20, 13.4 1, 1^.6 32.02,39.^2. 
38.22, 44.42, 50. 6J, 56.83. 63.04. 69.24, 

JLmjsiicai nL tmi = S.3MyT as a. fiuictboii oi Lin: viowiDg AOgt^i; 0v = 60° (top) and £)v = .lEQ" (bodtacn). 

Figure 7. Synthetic observations of the source with rj = 0.02 and Mach = 40. Left: polarisation maps. Vectors follow XB and their length 
is given by p. Vectors are superimposed on linear contours of // (/}. Right: logarithmic grayscale maps of // (/). 
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Figure 8. Synthetic observations of tiie source with ri = 0.02 and Mach = 80. Left: polarisation maps. Vectors follow XB and their length 
is given by p. Vectors are superimposed on linear contours of // (/). Right: logarithmic grayscale maps of // (/). 
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the magnetic energy in both the cocoon and beam over the 
source size for all runs. The curves are indeed strictly or- 
dered according to cocoon width. All the lighter jets have 
more magnetic energy than any of the light ones. Among jets 
with a given density, the faster ones have more magnetic en- 
ergy. Therefore, the underlying reason for the increase of 
the magnetic energy is the increase of the toroidal compo- 
nent due to the cocoon expansion, just as in Gaibler et al. 
(2009). There is another detail that confirms this finding: 
As described above, we find the usual 2D and 3D-phases for 
our jet simulations. The described amplification mechanism 
is very different in each phase. During the 2D-phase, field 
loops released in the jet head expand axisymmetrically, and 
substantial work is required to stretch them. Energy from 
this work is later found in the magnetic field. In contrast, 
during the 3D-phase the dentist's drill moves the jet head 
away from the axis in different directions. Field loops there- 
fore do not have to expand to reach large distances from the 
axis. They may keep their size, and get pushed into different 
corners of the cocoon at different times. Hence, once cocoon 
inflation reaches and goes thought the 3D-phase, almost no 
work is put into the fleld anymore. We believe this mech- 
anism causes the turnover in the magnetic energy seen in 
Figure O This turnover is visible for all the lighter jets at 
the comparable source size. The light jets, on the other hand, 
do not show much of an amplification in the first place, and 
also remain quite straight, i.e. essentially in the 2D-phase 
up to the end of the simulations. As expected, we do not 
find the turnover there. We note that a similar turnover is 
not found in the axisymmetric simulations by Gaibler et al. 
(2009) either, which is of course expected if it is linked to 
the 3D-phase. 

Why do we not see a linear increase with distance from 
the jet axis in the toroidal field like Gaibler et al. (2009)? 
Because of the 3D nature of the cocoon turbulence in our 
simulations. While axisymmetric turbulence can only stretch 
and compress a given toroidal field, 3D-turbulence may also 
turn toroidal field into poloidal one. The result is a turbulent 
cocoon field, with no geometrical similarity to the 2D-result. 

We see a strong axial field along the edge of the co- 
coons (see Figure [2|. This is due to velocity shear in this 
region (Section Q where the time average average back- 
flow speeds with respect to the ambient medium are about 
5xl0^ 5xl0^ 28xl0^ 22x10^ and 27x10^ km s-\ for the 
sources as they appear in Table 1. 

In the shocked ambient gas, on the other hand, mag- 
netic flelds are flrst compressed in the bow shock, and then 
reduced again due to adiabatic expansion of the gas, as it 
leaves the shock towards the cocoon. The effects of cocoon 
expansion on CMFs will be investigated in a sequel paper 
(Huarte-Espinosa, Krause, & Alexander 2011b, in prep.). 

The flow structure in our simulated radio sources is 
dominated by large scale motions, namely the toroidal and 
poloidal stretching mechanisms we have discussed in this 
paper. We cannot claim to represent the turbulence in our 
simulated cocoons well, because the resolution is too poor. 
Higher resolution should add additional small scale struc- 
ture, unless p revented by a sufficiently str ong magnetic field 
(compare e.g. I Krause fc Camenzindll200ll '). Yet, also turbu- 
lence is expected to have most power on large scales. There- 
fore, while higher resolution studies will still be useful, we 
would expect the results discussed in this paper to hold. 



3.2 Synthetic radio maps 

In Figures 21 [7] and [S] we present synthetic radio and po- 
larisation angle maps of four of our simulated FR II radio 
sources (all but the lighter-relativistic one) for three differ- 
ent snapshot times. These maps essentially reflect the field 
structure discussed above. The emission is dominated by fil- 
aments, hotspots and sometimes jets are seen. This is similar 
to what has been found in earlier studies, as detailed in the 
Introduction. The jet head region is more prominent at ear- 
lier times and for higher jet density. Our lighter jets feature 
more diffuse j et heads reminiscent of the sho ck web complex, 
described by iTregillis. Jones, fc RvrJ (|200ll ) . Our polarisa- 
tion angle maps are all dominated by larger patches. This 
is due to the fact that the cocoon dynamics is dominated 
by large vortices, about the cocoon radius in diameter. The 
backflow in our rj — 0.02 jets remains quite smooth and, con- 
sequently, the polarisation vectors are even more parallel to 
the jet axis than they are in the lighter -q sources. Generally, 
we find an almost one-to-one correspondence between the 
flow fleld and the polarisation vectors, as expected. 

In addition, in panel (b) of Figures IH (6] [7] and [8] we 
show synthetic radio and polarisation angle maps at differ- 
ent viewing angles. We consistently see that at small angles 
the axial fleld component gets smaller due to the projection 
effect, while other field components become relatively more 
prominent. 

The patchy distributions in our polarisation maps are 
in good agreement with typical observa tions of FR II 
sources and radio-loud quasars (see e.g. iBridle fc Perle\ 
1 19841 : ISaikia fc Saiteilll98a : [Gilbert et aLll2004l : lMullin et al 
200ii)- Along the projected direction of jets we see that 
|xb!<20°, which are smaller angles than elsewhere inside 
the cocoon. For rj =0.02, |xb| increases progressively along 
the vertical direction, from the jet axis to the edge of sources 
(Figures [7] and [8]). This is because in these simulations co- 
coons are narrow, and therefore the beam contributes signif- 
icantly to the emission, which is not the case in most of the 
observed sources. In contrast, for rj =0.004, |xb| shows weak 
trends along the vertical direction. The outermost vectors in 
all the maps are commonly tangent to the dimmest emission 
contours. This is similar to observations, but is of course in- 
fluenced by numerical problems at the contact surface, as 
outlined above. 

Polarisation degrees within 37-51 % are found inside co- 
coons, but higher, up to ~ 63%, both at the edge of sources 
and at the position of jets. We often see uniform patches 
with very similar values of |xb| and p at the position of 
bright emission shocks. The vectors frequently follow strong 
intensity gradients perpendicularly and have p > 50%. Re- 
gions of non-uniform |xb|, on the other hand, are fre- 
quently located between emission shocks (see e.g. Figure 2}. 
These co rrelations are in good agreement with observations 
fe.g. see lHogbomlll979l : lLaine"l9 8ll: iBridle fc Perlevlll984l: 
fSaikia fc Salter) ll988|r H a rdcastl e et al.l Il997l : iLeahv et all 
19971 : iGilbert et al.ll200'l_ and with models of p lasma com- 



pression and shear fe.g. ILaindlTgill : lMillej|l985l ). 

Our synthetic emission maps often show hotspots at the 
location of the jets' working surfaces as well as filaments in 
the radio lobes. Radio hotspots and filaments a re often seen 
in we ll resolved FR II sources (e.g. Cygnus A, iPerlev et al.l 
\1984 ). 
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We see the backflow of the (anti-paraUel) jets coUide 
and form sheets near the cocoon equatorial plane (the one 
normal to the jets and containing the central engine). There, 
our polarisation maps show B- vectors with |xb| > 7r/4 
above and below the centre of Figures |4}{8] (left column). 
We found instances of such polarisation angle distribu- 
tions in the observations of 3C 34 , 3C 336 and 3C 341 
dJohnson . Leahy, fc Garringtonlll995l . lMullin et al.ll200^ and 
[Gilbert e t al. 2004, respectively). 

At the end of the simulations we see laminar flows in 
the cocoons and also that both Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities are growing at the contact surface 
(Figure[T]). Such flows form tube-like structures or filaments, 
as we see in our synthetic emission maps (Figures |4]-??). 
Figure [9] shows emission maps of the right lobe of the model 
source with j) =0.004, Mach = 40 and tj^t =14.1Myr, for 
viewing angles of 30, 60 and 90 degrees. It is clear the struc- 
ture at the centre of the figures gets shorter and dimmer as 
the viewing angle decreases. This suggests a tube-like geom- 
etry for this feature. All our synthetic emission maps show 
filaments, however, we note that they form earlier in the 
low-r; sources than in the high-?; ones. 

3.3 Polarimetry and statistics. 

In order to analyse our synthetic observations we have pro- 
duced histograms of the polarisation angle and the degree 
of linear polarisation; see Figures [10] and llll respectively. 
In what follows we will see that these distributions show a 
clear correlation with the viewing angle, the jet-to-ambient 
density contrast and time too, but only a weak dependence 
on the jet velocity. 

3. 3. 1 The role of the viewing angle. 

The polarisation angle histograms are similar for all runs at 
0v =90°: They are all peaked towards zero degree, corre- 
sponding to the jet direction. The more isotropic distribu- 
tion at lower viewing angle is consistent with cocoon tur- 
bulence. We see only the distribution of the heavier jets 
remains peaked at a viewing angle of 60°, because of the 
weaker cocoon turbulence in these sources, relative to the 
ones with lighter jets. This confirms the magnetic field struc- 
ture is determined by the relative importance of turbulence 
as well as the amplification of the axial field due to the back- 
fiow in cocoons. 

For =30°, 60° and 90°, the mean value of |xb(»7 = 
0.004)1 is generally larger than that of |xb(7? =0.02)| (see 
Section [3.3.21 below). The dispersion of the polarisation an- 
gle seems to follow this trend as well. The differences are 
pronounced for viewing angles of 60 and 90 degrees and re- 
lated to the size of the data sample, i.e. the cocoons' vol- 
ume, which is inversely proportional to 77 in a non-linear 
way (Section l3.ip . Polarisation angle histograms at 6^ =30° 
show both the flattest gradients and the least number of 
vectors amongst all histograms, and their distribution does 
not show a Gaussian functional form. 

As the viewing angle increases we find the mean polar- 
isation angle, {|xb|), decreases non-linearly (see Figure fTU)) . 
(panel b, FiguresflHll). On average, IxbI diminishes for about 
9° for viewing angles from 30 to 60 degrees, and about 4° 



for viewing angles from 60 to 90 degrees. Cocoons have 
geometries that resemble prolate spheroids and thus mag- 
netic fields inside them should relax easier along the jet 
axis than towards the equator. However, to produce the syn- 
thetic maps we follow two steps: (i) rotate the sources anti- 
clockwise, perpendicularly to the jets, and (ii) project them 
onto the plane of the sky. Hence only the magnetic compo- 
nent along the jet axis (the horizontal one in the maps) is 
affected in this process and grows in proportion to cos(6v). 

The dependence of the degree of linear polarisation on 
the viewing angle is relatively modest and particularly evi- 
dent for r) =0.02 (Figure [TT|) . We see (p) increases about 7% 
from 6v = 30° to 9^ = 60° , and also about ~ 3% from 
= 60° to 6v = 90° . Hereafter, (p) represents the arith- 
metic mean of p. For all rj, the number of pixels in the po- 
larisation degree histograms consistently scales up with the 
viewing angle, in relation to the projected area of cocoons. 

3.3.2 The role of the density contrast. 

The jet-to-ambient density contrast is well known to be im- 
portant for t he evolution of cavities formed by astrophysical 
jets (see e.g. iKrausd l2003t ). Our synthetic maps show the 
density contrast also plays an important role on the radio 
source polarimetry. In general the projected area of sources 
is inversely proportional to r; in a non-linear way. Thus we 
see less magnetic fields in polarisation measurements with 
77 =0.02 than, 0.004. 

Given a timestep and a viewing angle, we find the mean 
polarisation angle is typically ~ 10° smaller for rj = 0.02 than 
for ri = 0.004. We see the spatial distribution of the polarisa- 
tion angle is more uniform for rj = 0.02 than for the lighter 
case. e.g. gradients greater than about 10° per computa- 
tional cell (10°/kpc) are less frequent in Figure [4] (left col- 
umn) than in Figure[71 corresponding to 77 =0.02 and 0.004, 
respectively. 

The statistical behaviour of the polarisation degree is 
very different. Given a timestep and a viewing angle, we 
frequently find higher values of the mean polarisation de- 
gree for T] — 0.02 than for the lighter case. On average, 
(p) {ri = 0.02)~ 47 %, while (p) (rj = 0.004)~ 42 %. Moreover, 
the polarisation degree histograms follow Gaussian-like dis- 
tributions. The mean polarisation degree at large viewing 
angles increases with time for the heavier jets, indicating 
that axial field line stretching gets even more important 
with time. Conversely, it decreases with time for the lighter 
jets, which shows that turbulence gets even more important 
with time for the lighter jets. The polarisation degree of the 
lighter jets does not depend on the viewing angle. 

3.3.3 Polarimetry evolution. 

The main features of the polarisation angle histograms seem 
to be shaped during the early expansion phase of the model 
sources, particularly for viewing angles of 60 and 90 de- 
grees. Here, magnetic fields tend to align with the jet axis 
(Xb =0°) as sources expand (panel a, FiguresflHHl). The con- 
sidered histograms decline steeply up to about 20°to 40°, 
and remain roughly constant for higher xb. The constant 
part is at a very similar level for all viewing angles of a 
given simulation. These findings correspond to the effects of 
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Figure 9. Synthetic emission filaments in the right lobe of the source with rj =0.004 and Mach = 40, at tj^t =14.1Myr. The structure 
at the centre of the figures gets shorter and dimmer as the viewing angle decreases from 90 to 30 degrees, from left to right, suggesting 
a tube-like geometry for the feature in question. 



isotropic turbulence, in combination with the stretching of 
field lines in cocoons, predominantly along the jet direction. 

The fractional polarisation evolves quite differently. At 
the first timestep (Figure [TT] black profiles), we see that the 
p histograms are fairly similar and show a linear relation, 
rising monotonically towards larger p. As the cocoons with 
r) = 0.004 grow, the profiles evolve into a peaked distribu- 
tion with a broad peak between about p = 0.2 and, 0.5. In 
contrast, the profiles for (77 =0.02, 0^ ^60°) always peak 
at p > 0.5. The fractional polarisation tells us about two 
things: (1) the degree of alignment of the field vectors that 
contribute to a given line of sight, and (2) the number of 
pixels long that line (assuming their contribution is differ- 
ent from each other). Hence we see the polarisation generally 
decreasing for 9 =30°. In cocoons where rj =0.004 we see 
a stronger p decline at early times than later on. This oc- 
curs because their expansion slows down at late times, as 
these sources approach pressure equilibrium with the ambi- 
ent medium. In cocoons where rj = 0.02, on the other hand, 
we see a slow sideways growth and thus p drops very slowly. 
Moreover, we see higher polarisation for larger ri, again re- 
fiecting that high density jets have more ordered magnetic 
fields and blow thinner cocoons. 

We note that an additional set of synthetic polarisa- 
tion maps (not shown) was produced assuming a spatially 
uniform distribution of synchrotron emitting electrons [i.e. 
Pc(x, t) = 1 in l|12p . for all x and t]. The polarimetric distri- 
bution of such maps was found to be very similar to the ones 
discussed in this section of the paper. Our results are not, 
therefore, sensitive to details of the electron distribution. 

Radio source polarimetry is related with the study of 
cluster magnetic fields because they induce Faraday rotation 
and depolarisatio n on the radio source emission (Pacholczvk 
1 19631 : lBurnlll966l ). Faraday rotation maps contain informa- 
tion about the ICM's magnetic structure (for a review see 
ICariUi fc Tavloij|2002t ). In a sequel paper (Huarte-Espinosa, 
Krause, & Alexander 2011b, in prep.) we will investigate the 
evolution of cluster magnetic fields using statistical analysis 
on synthetic RM observations which are produced using the 
expanding model sources we present here. 



4 DISCUSSION 

About a handful of studies on synthetic synchrotron emis- 
sion and polarimetry of extragalactic radio sources exist in 
the literature. We are not aware of any study that uses mag- 
netic fields evolved in a magnetohydrodynamic simulation 
self-consistently with the jet, as presented here. 

IJoned (|l988l ) modeled relativistic jets with a turbulent 
magnetic field ansatz, advected with the fiow velocity of the 
jet, to study the relation between linear and circular polar- 
isation in compact radio sources. The underlying hydrody- 
namic simulation is a conically expanding beam. With this 
ansatz, he gets a few, up to 23 per cent, linear fractional po- 
larisation. Though we start from a similarly turbulent field 
in our initial injection region, we get about 50 per cent lin- 
ear fractional polarisation in our beams. The reason is the 
order induced by the field line stretching described in more 
detail in Section 13.11 Models of jet coUimation and acceler- 
ation typically require a poloidal field near the source (e.g. 
IPorth fc Fendt , 2010. and references therein). The coherence 
length of this initial field should be small compared to ob- 
served jet sizes. Hence, stretching of the magnetic field in 
the beam seems to be an unavoidable consequence. The ef- 
fect is also found by Gaibler et al. (2009). The latter study 
is however the only one we are aware of that has employed 
a zero gradient boundary condition in the jet nozzle. 

Axisymmetric hydrodynamical simulations of coUi- 
ma ted light jets, similar to ou r approach, were employed 
by iMatthews fc Scheueij (Il990l ) to simulate the advection 
and deformation of passive magnetic fields s et up with a n 
initial isotropic random geometry, similar to IJonesI l|l988l ) . 
Matthews fc Scheuer implemented magnetic fields using pas- 
sive tracer particles and followed the distortion of the fiuid, 
at the respective position of tracer particles, by the veloc- 
ity field computed in the hydrodynamic simulation. This 
gives a reasonable approximation of the magnetic field struc- 
ture for dynamically passive magnetic fields. We also have 
a dynamically passive field, with a very similar initial, and 
slightly different nozzle boundary condition. In contrast to 
them, we do a full MHD treatment for the magnetic field. 
We confirm almost all of their results regarding the mag- 
netic field structure: Matthews fc Scheuer discuss in detail 
the toroidal and the poloidal stretching mechanism. As ar- 
gued above, we believe the toroidal stretching mechanism is 
mainly responsible for the magnetic energy increase in co- 
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Figure 10. Histograms of 1xb|- 9v increases from left to rigiit, column-wise. 



coons. We do not observe a dominant toroidal field directly 
because non-ajcisymmetric shear converts this component to 
poloidal field. We generally see a predominantly axial field 
component in cocoons, consistent with their poloidal stretch- 



ing mechanism. In contrast to them, we also find axially 
stretched and amplified magnetic fields in the beams. 

As iMatthews fc ScheueJ (|l990l ) do, we find field line 
stretching along the contact surface that separates cocoons 
from the ambient medium. In a resolution study they show 
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that the extent of that region gets smaller at higher reso- 
lution, but the field strength increases due to the increased 
shear. They also address synchrotron losses in the energy 
distribution of relativistic electrons. Due to such losses, they 
find that the aforementioned shear layer is very weak in syn- 
thetic radio maps. In our maps these features appear as edge 
enhancements and are likely a numerical artifact because our 
treatment does follow synchrotron losses. In reality, the two 
fluids may slip easily and the shear layer may be insignifi- 
cant. This depends on the magnetic viscosity of the plasma 
and is beyond the scope of this discussion. 

Moreover, as we do, iMatthews fc Scheued ljl99d ) find 
filaments in the synthetic emission images, but they report 
close to maximum fractional polarisation. We find lower, 
more realistic, fractional polarisation values in our simula- 
tions, especially for the lighter jets. The main reason for 
this difference is the breaking of axisymmetry: This allows 
for 3D turbulence in cocoons and for different directions of 
magnetic field vectors along the azimuthal direction. Yet, we 
also see fractional polarisation values in the cocoon body, 
far away from the beams and the edges, which are still 
somewhat high. This might be a resolution issue: the mag- 
netic field energy spectrum is close to Kolmogorov, which we 
have checked for the final snapshots of all our runs. There- 
fore, dominant structures are the large scale ones, which we 
should be able to capture. However, the roughly 50 cells 
we have over the fatter cocoons might still be too little to 
capture some important small scale structure that could re- 
duce the fractional polarisation. Our simulations show the 
fractional polari sation is very similar for di fferent jet ve- 
locities. Also, as iMatthews fc Scheuerl (jl990l ) have already 
noted, cocoon magnetic fields are largely independent of the 
initial conditions prescribed at the base of the beams. It 
is therefore unlikely that there is something fundamental 
to the cocoon structure that we miss. Another reason for 
low fractional polarisation might be that our jet densities 
may still be too high. We find the cocoon width, which is 
mainly regulated by the jet density, is an important factor 
for the polarimetry. Observed cocoons are usually wider rel- 
ative to the beam than the ones we produce here. This fact 
indicates lower jet densities in the observed radio sources 
l| Alexander fc Poolev.,1996. ; ^Krause 2003). Hence, low frac- 
tional polarisation might be yet another consequence of jets 
being very light com pared to their surroundings. Finally, 
IMatthews fc Scheueij (|l990) found small regions where field 
amplification and therefore synchrotron cooling became very 
significant in their simulations. In our 3DMHD simulations 
we see filaments in the cocoons (Figures 4-8, right column) 
and that magnetic fields there are about an order of mag- 
nitude str onger than the mean field. T hus we confirm the 
find ings of [Matthews fc Scheueij lfl990l ) . 

iTregillis. Jones, fc Rvul l|200ll ) carried out 3D-MHD 
simulations of a jet with r] = 0.01, Mach=80 and a helical 
magnetic field, the axial part of which extended throughout 
the computational domain. These authors studied the dif- 
fusive shock acceleration and transport of synchrotron rel- 
ativistic electrons. We d o not follow such processes. Then, 
iTregill is. Jones, fc Rvul l|2004 ) produced detailed synthetic 
observations of both the synchrotron and the X-ray -due 
Compton scatter from CMB photons. They emphasise that 
along the lines of sight that pass through strong shocks, 
most of the emission may come from regions close to the 



shock, and thus have close to maximum fractional polarisa- 
tion values. We might miss some of such regions due to the 
limitations of our simple model for the distribution of rel- 
ativistic electrons. However, the emission from the bulk of 
cocoons cannot be dominated by such features, as the frac- 
tional polarisation we predict for such regions is too high 
(compare above). This would mean that real radio lobes are 
relatively uniformly illuminated by relativistic electrons and 
are not dominated by relatively few isolated shock features. 



5 SUMMARY AND CONCLUSIONS 

We carried out 3D-MHD numerical simulations and syn- 
thetic observations to model magnetic fields in expanding 
FR II sources located at the core of a non-cool core galaxy 
cluster. A stratified fully ionized ICM was implemented, 
threaded by randomly tangled magnetic fields with a Kol- 
mogorov power spectrum. CoUimated, hypersonic and bipo- 
lar jets were injected in the centre of the computational do- 
main. The geometry of the jets' magnetic fields is initially 
random, and then shaped by the dynamics of jets. Jets form 
cocoons filled with light gas and magnetic fields, the struc- 
ture of which is determined by both the jets' backflow, via 
shear and compression, and the cocoon expansion. 

We have presented five simulations exploring the 
parameter space given by jet-to-ambient density con- 
trasts of rj ={0.004, 0.02}, and jet velocities of 
Wj ={40, 80, 130} Mach. We use the resulting model sources 
to produce synthetic synchrotron emission and linear polar- 
isation maps at viewing angles of 9^ ={30°, 60°, 90°}. The 
simulations have taught us the following. 

While we do not inject magnetic energy at the jet noz- 
zle, the magnetic energy in jets, and their host cocoons, in- 
creases with time. The amplification is stronger for wider 
cocoons, which are obtained for lighter and faster jets. 
The main amplification mechanism is the toroidal field line 
stretching ([Matthews fc Scheuedll990l : Gaibler et al. 2009). 
The toroidal field is however quickly converted to poloidal 
field and the resulting field structure is hence a competi- 
tion between MHD-turbulence and poloidal field stretching. 
Lighter jets are more turbulent and their magnetic field is 
therefore less aligned with the jet axis. 

Our synthetic polarisati on maps are in good agreement 
with radio observations (e.g. [ Johnson. Leahy, fc GarringtonI 
1 19951 : iGilbert et aLll2004i : iMullin et al.l 12000 ) We generallv 
see B-vectors that are parallel to the jet axis, tangent to 
the source boundaries and perpendicular to strong emission 
gradients. The degree of linear polarisation along both the 
jet axis and the source boundaries is higher than both inside 
and between radio lobes. 

The cocoon magnetic structure shows a strong relation 
with rj and a rather weak relation with Vj. In our polarisa- 
tion maps this occurs because the projected sources' area 
onto the plane of the sky is proportional to the cocoons' 
volume. The intrinsic polarisation angle distribution is con- 
sistently more uniform for rj =0.02 than in the lighter case. 
The mean polarisation angle is ~ 10° smaller when ij = 0.02 
than in the lighter case. Also, the intrinsic linear polarisation 
degree in the rj = 0.02 case is higher than in lighter sources, 
i.e. when rj = 0.02 we see p within 46-51 per cent in the co- 
coons and, ~ 63 per cent at the sources' edges. Conversely, 
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when rj — 0.004 we see p within 25-45 per cent in the co- 
coons and, ~ 63 per cent at the edges. Even for our hghter 
cocoons, the fractional polarisation is somewhat high away 
from the edges and beams, which might be a resolution issue 
or due to the fact that our cocoons are thinner than those of 
most observed FR II radio sources, which is related to the 
jet density. 

The distribution of the polarisation angle (magnetic 
vectors) depends on the viewing angle between jets and the 
line of sight, 6^. On average we see (|xb|) decreases about 
9 degrees as 9^ goes from 30° to 60°, and about 4 degrees 
as 6v goes from 60° to 90°. In contrast, only {p (rj =0.02) ) 
shows an increase of about 7% as 9v goes from 30° to 60°, 
and also about 3% as 6^ goes from 60° to 90°. This is be- 
cause cocoons have geometries similar to prolate spheroids, 
inside which the poloidal momentum flux is higher than 
the toroidal one. Cocoon magnetic fields are thus mainly 
stretched along the polar direction (the jet axis) which pro- 
jection onto the line of sight is proportional to cos(Sv). 

We see the main features of the |xb| histograms are 
shaped during the early expansion phase of sources, partic- 
ularly for 6v > 60°. In this case, magnetic fields tend to align 
with the jet axis as sources grow. For 6^ — 30°, on the other 
hand, xb is distributed nearly isotropically. The fractional 
polarisation is broadly distributed around about 30-40 per 
cent, and decreases in time. 
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